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Abstract 

In the framework of irreversible thermodynamics, we studied the transport prop- 
erties of QGP. Shear viscosity and non-equilibrium entropy density related to viscous 
process at finite density has been investigated in weakly coupled limit by using ki- 
netic theory. The results show that the chemical potential increases viscosity yet 
decreases the non-equilibrium entropy density and thus contributes positively to 
their ratio compared to the pure temperature case. As to the temperature depen- 
dence of the ratio, the ratio first decreases rapidly and then increases in the physical 
region, presenting a minimum value of 0.4 at the temperature around 182MeV. 

1 Introduction 

Recently scientists believe that the quark-gluon plasma(QGP) found at relativistic heavy 

ion collider(RHIC) is strongly coupled, contrasting with the weakly coupled as expected 

before, which is the so-called sQGP [1-3]. They also believe when the temperature goes 

higher, for example T > 2 — 4T C , the hot medium would be expected to become a weakly 

coupled system with some dissipative structures. Additionally, investigations on compact 

star demonstrate that the viscous properties are significant in explaining many of their 
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behaviors [4,5]. Therefore dissipative structure of strong interaction system, especially 
the sQGP, is a remarkable topic. 

In irreversible thermodynamics, the discussions on the dissipative properties of a sys- 
tem are focused on the entropy production in an unit time 



where Xi is the thermal force which is determined by the gradients of energy, temperature, 
chemical potential etc. 7$ is the corresponding flow driven by X,i which can be written in 
the linear response approximation as 



where Ly are the transport coefficients. Thus it can be seen clearly that the entropy 
production is determined by two factors: one is the thermal force Xi which is the the 
external cause describing the environment; the other factor is the transport coefficients 
Lij which are the intrinsic causes reflecting the responsibility of the system driven by 
the thermal force. Generally speaking, the entropy variation AS in an unit time can 
reflect the evolution of the dissipative non-equilibrium state. Therefore the Eqs. (JTJ) 
and (J2J) provide us the basic evolutional information of a dissipative system theoretically. 
Inserting Eq. (J2J) into Eq. ((TJ, one could also analyze the state, or namely the thermal 
force, of a dissipative system through the ratio of transport coefficients to the entropy 
production which might be inversely obtained through fitting the experimental data. As 
for the QGP produced at RHIC, this general discussion is focalized on the shear viscosity 
to entropy density ratio which determines whether it is feasible to describe the hot and 
dense medium by the ideal energy- momentum tensor Tq V . Besides the shear and bulk 
viscosity induced by the velocity gradient, this tensor is also relevant to the thermal 
conductivity induced by the energy gradient, which is zero with Landau and Lifshitz's 
definition on the hydrodynamic velocity [6]. 

In heavy ion physics, the evolution of the QGP density is almost longitudinal, with ve- 
locity gradient along the transverse direction. Thereby bulk viscosity is usually neglected 
due to the domination of shear viscosity. More proof from recent lattice simulation favors 
the ignorance of the bulk viscosity since it is much smaller than the shear viscous one [7]. 
As a result, the shear viscosity and entropy density ratio becomes a very important means 
of understanding the dissipative properties of the QGP produced in heavy ion collision. 
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As far as the sQGP is concerned, Ref. [7] demonstrated a result from lattice Monte 
Carlo simulation, pointing out the ratio of the shear viscosity to the entropy density is 
smaller than unit but most probably larger than the universal bound 1/An which is ob- 
tained from the gauge theory /gravity duality [8-11] in equilibriate superstring theory. 
While one must notice the entropy density in this ratio is supposed to be an equilibrium 
one, which is not exactly the one appears in the Navier-Stocks equation. Some other 
estimations [12], e.g. from the Heisenberg uncertainty principle, has given a result as yg. 
Studies on this strong coupled mechanism are going on since it is desirable to understand 
the perfect behavior of the production at RHIC. As to wQGP (weakly coupled QGP), 
where the perturbative expansion works, people knew much more about its viscous prop- 
erties than those of the sQGP. The transport coefficients, especially the shear viscosity, 
were discussed by many authors in weakly coupled limit in the kinetics theory [13-21] and 
in the thermal field theory via the Kubo formulae [22-33]. From the theoretical points 
of view based on the previous analysis, the ratio of shear viscosity to entropy density is 
also the foundation of understanding the wQGP's dissipative properties. However most 
of those investigations mentioned above are usually focused on the extremely high tem- 
perature but zero chemical potential environment. Further studies on the ratio of shear 
viscosity to entropy density at finite density should be carried on, which is one of our 
motivations. In addition, as for the denomination of the ratio, only the equilibrium en- 
tropy density [8-11,34] and the entropy production in phase transition [12] are discussed. 
The entropy production induced by the dissipative forces should be included. In this pa- 
per, considering these two factors we calculate the shear viscosity of the wQGP at finite 
chemical potential and the corresponding entropy production in dissipative processes self- 
consistently so as to obtain the viscosity to entropy density ratio, representing a upper 
bound in the weak coupling limit. 

The paper is arranged as following: shear viscosity and entropy density are evaluated 
in the framework of kinetic theory in Sec. 2 and 3 respectively. Summary and discussions 
will be presented in Sec. 4. 

2 Shear viscosity 

In this section, we derive an expression for shear viscosity rj of hot QGP at finite chem- 
ical potential in leading logarithm order in the framework of kinetics theory by using 
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variational approach, following Arnold et al [19] and our previous paper [21]. 
2.1 Formalism and definition 

In a state that is departed not far from the equilibrium, the energy momentum tensor 
can be decomposed into ideal and dissipative parts as 

= Tg u + tt"" = (e + P)u»u v - Pg^ + tt"" (3) 

where e, P, are the energy density, pressure, and viscous shear stress, respectively. 
The four-velocity u^(x) in the local rest frame is (1, 0, 0, 0). In the linear response theory 
with Landau-Lifshitz convention, the fluctuation of spacial stress tensor is proportional to 
the first order of velocity gradients with neglected bulk viscosity by defining the coefficient 
i] as shear viscosity, 

8 fa) = -v(9iUj + djUi - \5ijd k u k ) = -rjXij, (4) 

where X^- is the so-called driving force. 

On one side, in kinetic theory the energy momentum tensor is the secondary moment 
of the distribution function, 

(Tij) = I ^ [g f f(t, x; p) + g f f(t, x; p) + g b b(t, x; p)] , (5) 
J p P 

where the momentum space integration j p is a shorthand for J A^s, and f(x; p), f(x; p) 
and b(x; p) are one particle fermion, anti-fermion and boson distribution function respec- 
tively, with gf,gj and g b of their degeneration degrees. 

On the other side, when one decomposed the one particle distribution function as 
a local equilibrium part plus a fluctuant part, namely f s = n s + 5f s where n s is the 
equilibrium distribution function and s denotes for species of particles, the fluctuation of 
the energy moment tensor, which is exactly contributed by the fluctuation of the stress 
tensor, is 

8{Tij) = [ — (g f Sf + g f Sf + g b Sb) = 6 fa) = -v^u, + d jUi - l5 t3 d k u k ). (6) 
Jp Po 

In principle, one could evaluate the shear viscous coefficient as long as the distribution 
functions are known. While it is not an easy task because the one particle distribution 
function satisfies the Boltzmann equation of the usual form 

(d d d \ 

- + p • — + F ext • - f s (t, x; p) = -(Cf s )(t, x; p), (7) 
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where F ex t is the external field and the C is the collision operator. Even in a simple 
steady system without external field, this evolution equation contains both differential 
and integral terms, where the collision integral might be rather complicated in different 
processes. In our paper, for a upper limit estimation, we just consider 2 — > 2 elastic 
collisions to obtain the leading logarithm result, although for the relaxation and the shear 
viscosity, the inelastic scattering might be significant especially in the initial stage of QGP 
formation. Under this assumption, the collision operator is defined as 

< c « p> - 5 L, IZ^m mp +K - pl - K,) 

x {f s (p)f s (k)[l±f s (p')][l±f s (k')} -/ s (p')/ s (k')[l±/ s (p)][l±/ s (k)]}, 

, where p, k, p' and k' denote the momenta of the incoming and outgoing particles re- 
spectively. \M. | 2 is the two-body scattering amplitude. The 1 ±/ s factor is the final state 
statistical weight, for boson with the upper sign and for fermion with the down sign. 
For convenience, one can rewrite the fluctuation of distribution function in the form 

Sf 8 = n a (p)[l ±n s (p)}ip s (x;p) (8) 

Inserting the decomposed distribution function into the right hand side of the Boltz- 
mann equation ((7j), one notices that (Cn s )(t; x, p) = when n s takes the form of Jiittner 
distribution as 

^(t;x,p) = e/3(p , u _ Ats)±i (9) 

with fx s = ±/i for both fermion and anti-fermion. Notice that in the local rest frame, 
this one particle distribution function is degenerated into the ordinary fermion and boson 
distributions, 

= e P(P-l) + V n ^ p) = e^) + V nb{p) = ^~l (10) 
Linearizing the collision term, one can obtain 

(C/ )(x;p) = \f ^|M a t(P,ir;P',^)rx(27r)V 4 )(P + ^-P'-^') 
1 J p'W bed 

x n a {p)n h {k)[l±n c (p')\[l±n d {k')\ 

x [^o(x; p) + (f b (x; k) - v? c (x; p') - y^(x; k')] , 

(11) 
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where abed represent the species of the particles and \M%~(P, K; P', K')\ denotes for 
(2pi")( v I )(2koY2k' ) ■ ^ ne sum m f ron t °f the matrix element means all possible collision pro- 
cesses relevant to the leading-log contribution are involved and properly treated without 
double counting or multi-counting. 

As for the left hand side of the Boltzmann equation, the gradients acting on Sf s give 
higher order in the departure from equilibrium, so that in the first order approximation, 
only n s should be considered on this side, namely, 

LHS = [3n s { P )[l ± n a (p)]/ -(p)^-(x) (12) 

with iy(p) = \ijPiVj — where the time derivative and the the external field terms 

vanish for viscosity and the spatial tensor X*j(x) is just the driving force defined in Eq. 
©• 

Comparing both sides of the Boltzmann equation, one finds that y?(x; p) on the right 
hand side must have the same angular dependence as the driving term, i.e. 

y? s (x; p) = (3 2 I ij (p)X ij (x.)xs(p) (13) 

where Xs{p) is rotationally invariant function depending only on the amplitude of mo- 
mentum. 

With equations (fTTj) . (fT2" j) and (|T3j). one cancels the driving field on both sides and 
recasts the Boltzmann equation into 



where 



s*.( P ) = (c*y( P ) (14) 

S£(p) = -T P n a (j>)[l±n,(jp)]Iii(p), ( 15 ) 

Xtj(p) = hj(p)Xs(p)- (16) 

To solve this equation, we introduced the variational approach following the basic 
steps in reference [19]. First, one defines an inner product as 

(f,g)={3 3 [ /(pMp). (17) 
Jp 

Then taking product on both sides of the Boltzmann Eq. (JT3)l and defining the functional 

Qbc] = (xij, Sij) - lixij, Cij) (is) 
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Figure 1: The possible processes which contribute to the leading- log in the collision 
term in QCD plasma. The solid line is for quark and the wiggly line is for gluon. 
The matrix elements arise from all the five diagrams are \Qg i dAC\ (3 — p — ^| — , 
8g*d F C F C A (^), 8g^ + 8g±d F C F (± + f) and -8^Cf + f) 

respectively, with ci^ = 8, dp = Ca = 3 and Cp = | for SU(3) group. 

with explicit form of each term as 

(Xij, Su) = -p 2 J2 I P /o(p)[1 ± IS (p)]X°(p). (19) 

a J P 

(XV,CXV) = ^/ ^|M^| 2 (27r) 4 5( 4 )(P + ir-P'-^) 
» ^P,k,k,k' afccd 

x n (p)n 6 (k) [1 ± n c (p')] [1 ± n d (k')] 

x [x&CpJ + xtW-^pO-xSCJO] 2 . ( 2 °) 

one may relate the maximal value of Q with the viscous coefficient through Eqs. (jUJ), (jEJ) 
and (|13|) with the fact that the maximum value of Q[x] occurs when x(p) satisfies Eq. 
(Hl,i.e., 

V = TsQmax\ x=Xmax - (21) 

2.2 Collision terms 

The collision terms of the Boltzmann equation are associated with different interaction 
processes. In the leading logarithm order, only five diagrams contribute, which are pre- 
sented in Fig.((TJ). 
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Before starting calculation of these collision terms , some arguments and kinematic 
relations should be manifested first. 

• Approximations 

Two approximations have been employed in our further evaluation. First, the high 
temperature approximation. We assume the temperature T of the QGP is extremely 
high which is the only 'hard scale' in the system, with other quantities like chem- 
ical potential /i, thermal mass etc. indicated by gT, are much less than the hard 
scale in weakly coupled theory. With this high temperature approximation, the 
fermion(anti-fermion) distribution function can be expanded in terms of fi/T, ne- 
glecting the thermal mass term directly [21]. The second important approximation 
is the forward scattering approximation, i.e. the momentum transfer between the 
incident particles q ~ gT is rather small so that the momentum difference between 
the incoming and outgoing particles on the same interaction vertex can be ignored 
in the distribution functions. 

• Kinematics 

We mark the momenta of the particles as shown in Fig^c). The three usual Man- 
delstam variables are defined as s = (P + K) 2 , t = (P — P') 2 and u = (P — K') 2 . 
With the forward scattering approximation, the incident angle between p and k 
can be described by 9 and 0, which are the angle between p and q and the angle 
between the pq plane and p'q plane, 

co$9 pk = 1 + (1 -cos 2 #)(l -cos0). (22) 

Now let us turn back to the Eq. ([20)1 . Performing the integral over dk' with the help 
of S 3 (p + k + p' + k'), one may transform the equation into 



'3 roc pq poo poo p2tt 



i.Xij ) CXij) 




\M%\ 2 n a (p)n b (k)[l ± n c (p))[l ± n d (k)) 

abed 

x [xU?) + XUV ~ *W) ~ ■ (23) 

where u is a dummy integration variable first introduced by Baym [17] with p' = p + uj 
and k' — k — u. 
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It is easy to discover that the integrand of Eq. (J2~3j) is composed by three parts: the 
matrix element, the distribution functions and the x term. According to the exchanged 
particle, the \ term can be sorted into two classes: Fig^a-c) belong to one category 
where the on-shell particles interact with each other through exchanging a boson, namely, 
they have the same species of incoming and outgoing particles on one interaction vertex; 
Fig^d) and (e) belong to another category where they contain off-shell fermions, i.e., 
different species of on-shell particles are bounded to the same interaction vertex. For the 
first category, the x term contributes a small q 2 in the forwarding scattering approxima- 
tion, which softens the infrared singularityand gives the leading-log form of viscosity [19]. 
For the two on-shell line of quarks, the \ term thus is specified as 

[XW) ~ 4(P)] 2 « [q " V4(P)] 2 = ^V(P)'] 2 + 3^^[X 9 (P)] 2 . (24) 

where x 9 {p)' = dx 9 {p)/dp. We assumed that the quark and anti-quark have the same 
departure behavior from the equilibrium state which is denoted by x q i since nothing 
but the viscous process is involved. The same trick will be performed on x 9 ■, the gluon 
X function, and leads to similar result. One can prove that the interference term like 
[x 9 (p) — X 9 (p')] [x 9 (k) — x 9 (k')] vanishes when carrying out the du and dip integration. 
Following the above discussion, we can list the all the x terms for different diagrams in 
Tablf^ where the only survived formats are (x q — X q ) 2 1 {x q ~ X 9 ) 2 an d {x 9 ~~ X 9 ) 2 - 

There are two facts that should be clarified. One is that all possible channels for 
one diagram are not to be added directly. In fact for each channel, one must perform 
a convolution integral for the exact x term from the very channel with the correspond- 
ing distribution function. The other fact to be noticed is the relationship between the 
Mandelstam variables under the forward scattering approximation. For instance, as to 
the massless on-shell particles, s = (P + K) 2 = 2P ■ K « — (P — K') 2 = —u in u- 
channel. Similarly, one has s ~ — t in t-channel. With this convenience, combining with 
Eq. (}2*2*|) one finds out that all the matrix elements reduce to only two forms, 

V + ^ f S > + u*\ -.*"*",! _<*»#)', (25) 



u 2 \ t 2 a A 

/ u— channel V / t— channel " 

t\ fu\ 2pk 



,2 



l-cos0). (26) 



. U J u-channel channel q 

The last part of collision integrand is the distribution function, which is the main 
distinguished feature for each diagram and channel. Analyzing each possible process 
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Processes 


X functions 


FigUTa) 


[x 9 (p)-X fl (p')] 2 + [x 9 (k)-x 9 (k')] 2 


FigCtbl) 
FigHi;b2) 


[x 9 (p)-x 9 (p')? + ix q (V-x q (V)] 2 

[x q (p) - X q {p')f + [X s (k) - X *(k')] 2 


FigHTc) 


[x 9 (p) - X 9 (P')] 2 + [X 9 (k) - X 9 (k')] 2 


FigtUdXe) 


[x ? (p) - X 9 (P')} 2 + [X 9 (k) - X 9 (k')] 2 



Table 1: The x terms for five diagrams. The Fig^b) has two sets of x functions be cause 
it involves different channels which bring on different momentum dependence of x q an d 
X 9 - 
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Processes 


Distribution functions 


Figffla) 


2n b (p)n b (k) [1 + n b (p)} [1 + n b (k)} 


Figntbi) 
Figiu;b2) 


2N f {n f (k)[l -n f (k)) +n f -(k)[l - n f -(k)]}n b (p)[l + n b (p)} 
2N f n b (k)[l + n b (k)]{n f (p)[l - n f (p)} + n f (p)[l - n f -(p)]} 


Figure) 


2N}{n f (p)n f (k)[l-n f (p)}[l-n f (k)]+n f (p)n f (k)^ 

+n f (p)n f (k)[l - n/(p)][l - nj{k)\ + n f -(p)n f -(k)[l - n/(p)][l - n/(fc)]} 


FigHId+e) 


2iV / {n / (p)n 6 (A;)[l - n/(A;)][l +n 6 (p)] +n / (p)n 6 (A;)[l - + n 6 (p)] 
+n&(p)n/(/c)[l - n/(p)][l + n 6 (A;)] + n 6 (p)n/(A;)[l - n/ (fc) ][l + n b (k))} 



Table 2: The distribution function terms for the five processes, where Nf is the quark 
flavor. The factors in front of the distribution functions are the freedom of degeneration, 
which are related to the distinguished reaction channels. For example, qq «-> qq appears 
4Nf times in the sum over species. 

carefully, one can obtain the distribution function term for each diagram in Fig. ^ which 
is listed in Tabel |2] 

Inserting the expressions in the three tables into Eq. ()23|) and carrying on the inte- 
gration over dk, d(f> and dou, one can finally obtain the full collision term, 

iXijt^Xij) = {Xiji^Xij) ^ {Xiji^XijY ^ "I" {XijjCXij)^ ^ 

+ (x lv C Xlj ) ic) + (Xv,C^) id+e) (27) 

where 

n 2j /T2 rT Art r°° 

(W.Crf''^ / - / dpn b (p)[l + n b (p)]{p 2 i X °(p)'] 2 + Q[x 9 ( P )] 2 }, 

671 Ja s T 1 JO 

(28) 
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{Xi . Cxi . r) = 2a,d F N l c F c A r * r» ^ {p2 W{j>)f + 6 m]2 

611 JolsT Q. JO 

x {n f (p) [1 - n f (p)} + n f {p) [l - n f (p)} }, (29) 

x {p 2 [x 9 (p)'] 2 + Q[x g (p)] 2 }, (30) 

x { n/ (p) [1 - n f (p)} + n f (p) [1 - n f (p)} }, (31) 

(x,,Cx^ = 2a2 fff^ f * rdpp^ip)-^)? 

^d A -K 6 J asT q J 

x {n,(p)[l+n b (p)] (y- 0-616^ + |^ 
+n f (p)[l + n b (p)} (y + 0.616| + ^ 
+b f (p)[l-n,(p)] (y + 0.616| + ^) 
+ n 6 (p) [1 - n f (p)] - 0.616^ + ^) } , (32) 

with a s as the fine structure constant. We replace the limits of dq integration by the 
hard and the soft scale T and a s T due to the leading- log treatment, which includes the 
self-energy effect of the exchange line [19]. 

2.3 Variational method 

So far we have obtained the right hand side of the Boltzmann equation, and the other 
side can be written down directly as 

o2 rcc 

(Xij, Sij) = 2 / dpp 3 \ d F N f {n f (p)[l - n f (p)} + n f (p)[l - n f (p)}} X q (p) 

71 Jo K 

+d A n b (p)[l + n b (p)] X 9 (p)} (33) 

In this subsection, we expand the trial function xip) by a finite set of basis with 
independent variational parameters to maximize the functional Q[x\- 
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As we have known that in the case of viscosity, the % function has two components, 
X q and x 9 , 

^ - ( S£ ) • 

Expanding the two components by the same basis 

N N 

X 9 (P) = ^2 a ^rn(p), X Q (P) = ^2 a N+m<Pm{.P), (35) 
m=l m=l 

one could read out the basis-set components of S and C 

ifiij; Xij) ^ Q'm'S'm; (X?J? ^Xij) ^ ^ C mn Q, n . (36) 

When the Boltzmann equation is satisfied, the coefficients in front of the bases will be 
expressed as a = C^S. Accordingly, the shear viscosity becomes 

^ = J^Qmax = if; - ' S = j^S C S (37) 

With the natural one function ansatz <pi{p) = p 2 /T, one can evaluate the integral of 
Eqs. (|27jl and analytically, 



7T 2 



(3 3 f°° d A n b (p)[l + n b { 

S = -—I dpp | 

d F N f {n/(p)[l - n f {p)\ + n f (p)[l - n f {p)}} 



120C(5)T 3 f d A \ 



n 2 



(a 



167r 3 «2T 3 In a," 1 
9^ 



' d A C A \d A C A 

+N f d F C F (l + ^)] 



\N f d F C F [d A C A {\ 

+ Nf d F C F (l + ^^)} J 



30 ji ■ 
7vr 2 T 2 - 

,2 



+ 



97T 2 

128 



N f d F C 2 F d A [ 1 - 0.013 




(39) 
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where ((x) is the Riemann ^-function. Then one can compute the coefficient vector and 
the shear viscosity of 2-flavor QGP at finite temperature and chemical potential as, 

, a u 1 ( -3.281 + 0.940^ \ , , 

V \ -6.792 + 0.199^ j 

and 

1.09T 3 _„/z 2 




a'ln^V™^ <41) 
which recovers the result in Ref. [19] when the chemical potential vanishes. 



3 Entropy production in viscous process 

From the non-equilibrium thermodynamical points of view, the entropy S symbols a scale 
that how far the system is apart from the equilibrium state where the entropy takes 
the maximum value. Other than the viscosity which reflects the intrinsic property of 
a system, the entropy depends further on the environment the system resides in, or the 
velocity gradient in our case. It is the velocity gradient that determined the evolution from 
non-equilibrium to equilibrium in a given system. Unfortunately, we know little about 
the velocity distribution which depends on the experimental environment and should be 
obtained by fitting data. In Sec. 4, a maximum estimation of the velocity gradient will be 
employed so as to demonstrate a maximum entropy production which leads to an upper 
bound of viscosity to entropy density ratio. 

We introduce two different schemes to calculate the entropy production and find out 
that their results are comparable qualitatively. 



3.1 Scheme I 

In relativistic kinetic theory, the entropy density of a system is defined as 

* = -Z)/ ^{A(x;p)ln/ S (x;p)+[1±/ S (x;p)]ln[l±/ S (x;p)]} (42) 

where g s is the degeneration degree of each particle species. 

Decomposing the distribution functions and inserting Eq. (jSJ into Eq. (|42j) . one finds 
that only the second order in ip survives while the first order in ip s vanishes due to the 
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rotational invariance of n s (p). The entropy density thus becomes 

s = s ° ~ / Yl 9sTls (p}[ 1± ns (p)] ( x ; p) ( 43 ) 

where s is the entropy density in the equilibrium state, 

s = -J {g f [n f lnn f + {l - n/) ln(l - n f )] +g f [n f \nnj 

+ (1 - nf) ln(l - nf] + g b [n b \nn b - (1 + n b ) ln(l + n b )] j 
= 16.22T 3 ^1 + 0.123^ (44) 

with gf — g^ — 2Nfdp and = 2<i^ and the integral evaluated by expanding the fermion 
distribution functions in terms of small [i/T. It is clear that the second term of Eq. 
()43|) denotes for the departure of the entropy density from its maximum value, i.e., the 
entropy production. As general, the non-equilibrium entropy density has the structure of 
s = s — As, as we expect. 

In this scheme we are going to directly evaluate the second term of Eq. (|43|) by 
adopting the solutions of the Boltzmann equation in Eq. ()40|) . leaving only the velocity 
gradient to be estimated. Inserting expression (fTB^l into Eq. (|4^j) with the solved x{p) m 
the variational approach, we can obtain 

/3 6 f°° 

s = s - j—^ dp p 6 {g f [n f (l - n f ) + n f (l - nf)] a\ 

= -J&f^O^U^V. 1 (46) 



(ajUna,, 1 ) 2 \ T 2 J \ dx 

It's obvious that the second term of Eq. (|43jl is the entropy production relative to 
viscous process when the system evolving to its equilibrium state, which depends on the 
kinetic parameter a s , the thermodynamic parameters T \i and the transport environment 
du z /dx. To understand further on their physical meaning, we introduce a second scheme 
to evaluate the entropy produced in the viscous process. 



15 



3.2 Scheme II 

From Scheme I, one can see the non-equilibrium entropy density is obtained by subtracting 
the entropy production from the equilibrium entropy density, which is denoted by 

s = s - J dXd^ = s - J dtd t s° (47) 

where s M = su^ 1 is the four entropy flow. The last equality is based on the description of 
longitudinal expansion model in local rest frame. 

As we stated in Sec. 1, the entropy production could be obtained by the product 
of the driving force and corresponding thermodynamic flux. Generally speaking, the 
entropy produced in unit phase space <9 M s M takes a bilinear form: it is the sum of all 
kinds of transport flows, each multiplied by a characteristic thermodynamic driving force. 
Furthermore, under the linear law, the transport flow is proportional to the the driving 
force with the coefficient of viscosity, conductivity, etc. Therefore, the universal expression 
for the entropy production should be the combination of squared driving forces coming 
from various transport processes, weighted by the corresponding transport coefficients. 
In our case, we would like to know the entropy production just from the viscous process, 
which is [13] 

.o _ V V 2 _ 2 V ( du z 



= d t s u = ±Xi» = 9f ^ (48) 

where is defined in expression Noticing that the integration variable t = rcoshy, 
where r describes the time scale of system departure, one could replace it by the relaxation 
time t v for the maximum estimation. With this approximation, one can perform the 
integration on the pseudorapidity plateau, 



TR\h 



cosh j/o 



^ so-^y-/ (sinh^coshy) (49) 



1.875T / u 2 \ (du z . 



a^lna^, 1 ) 2 \ T 2 J \dx 

where y = 2 is the edge value of the pseudorapidity plateau and will be calculated in 
the appendix. 

We argue that this estimation is consistent, for the relaxation time defined in the 
Boltzmann equation describes the one particle property in the proper reference frame. 
From Eq. ()48|) one can see clearly that the entropy produced in viscous process depends 
not only on the intrinsic property r\ but also on the exterior transport environment of 
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the system. These internal and external factors determine the dependence of the entropy 
density on the kinetic and thermodynamic parameters, and also the velocity gradient. 

4 The ratio of 77 to 5 

Now we are going to turn to a significant quantity i]/s which has been mentioned in 
the first section. In a pure viscous transport process, the state deviating from the equi- 
librium is solely stimulated by the spacial derivative of velocity field d{Uj. For a rough 
estimation of the velocity gradient for a longitudinal Lorentz boost invariant system, 
we adopt Brioken's assumption = x^/r [35], where x^ is the space-time point and 
r is the proper time. In this scenario, the shear velocity is just along longitudinal di- 
rection z which is perpendicular to the velocity gradient direction x. Thereby diUj 
has just one non-zero component du z /dx. Noticing that in the light-cone coordinate, 
u z = x z /t = sinhy where y denotes the pseudo-rapidity and the maximum of velocity 
gradient is du z /dx ^ u z /Ra = sinhy / Ra when y takes the edge value of the pseudo- 
rapidity plateau. Here Ra is the nucleus radius and the system is considered rotational 
invariant around the z-axis. In Au+Au 200GeV collision at RHIC, the pseudo-rapidity 
plateau is from — 2 ~ 2 [36] in the most central collision, and the radius of gold nucleus 
is about Zh.Z&lGeV- 1 . 

With this maximal velocity gradient estimation, the ratio of viscosity to entropy den- 
sity obtained in the two schemes are 



(52) 

respectively. These two results are qualitatively comparable for they have exactly the 
same dependence on dynamic and thermodynamic parameters except some differences on 
the factors. 

We manage to demonstrate the basic features of the dependence of the ratio on the 
thermal parameters by adopting the result from scheme I. As their chemical potential 
dependence is quite clear, we here just present their variation on the temperature. To 



(51) 



a) > a 

^ so 



0.23GeV 2 



1 + 0.12^ 
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Figure 2: The shear viscosity to entropy density ratio of 2- flavor QGP. 



isolate the temperature dependence, we use a simple form of the running coupling constant 



where T c is the critical temperature of QGP phase transition. We assume the chemical 
potential has negligible effect on the constant since it is small compared to the temperature 
above T c . Inserting Eq. (p)3*j) into Eq.([51|). we can obtain Fig. 2 with T c = 175MeV and 
/I = 46MeV. 

It should be noticed here that in Fig. 2 the temperature dependence of the ratio 
has been extrapolated into the strongly coupled region where a s > 1 to see where is the 
physical boundary and the whole variety of the ratio in the physical region. One may 
find that the structure of the viscosity to entropy ratio is rather complicated due to the 
appearance of two singularities, which are contributed by the logarithm and the denom- 
inator of the enhancement factor. To distinguish the physical curve from the complex 
structure, we naturally introduce a physical criterion: both the ratio in equilibrium and 
the enhancement factor, which are presented in Figs. 2 and 3, should be positive. 

It is clear in Figs. 2 and 3 that the physical curve is the one on the right side of the right 
singularity, i.e., the physical region is bounded by a temperature, which varies according 
to the constant over the fine structure constant under the logarithm. With detailed 
observation on the physical curve, one finds that it does not decrease monotonically 
with the temperature but presents a minimum value of about 0.4 at T = 182.1MeV^ 
instead, which is twice greater than the phenomenological estimation [37]. This minimum 
value, seeing from mathematical point of view, is owing to the competition between the 



[7] 




(53) 
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Figure 3: Ratio of shear viscosity to entropy density in equilibrium. 




Figure 4: Enhancement factor. 
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weakening running coupling constant and the rising temperature itself. 

5 Summary 

In the framework of irreversible thermodynamics, the transport properties of so-called 
QGP produced in heavy ion collision have been studied. Based on the Boltzmann equa- 
tion, we calculated the shear viscosity of two-flavor QGP at high temperature and finite 
density in weakly coupled limit in kinetic theory. The result shows the finite density effect 
provides positive contribution to the shear viscosity by adding a quadratic term of \ijT to 
the pure temperature case. Furthermore, we calculated the non-equilibrium entropy den- 
sity at finite chemical potential through two different schemes by subtracting the entropy 
production due to shear viscous process from the final equilibrium entropy density. In our 
evaluation, the thermal force in the shear viscous process, namely, the velocity gradient, 
was estimated by it maximum value considering the longitudinal expansion in heavy ion 
collision. Finally the viscosity to entropy density ratio was demonstrated in a extensive 
temperature regime including both physical and unphysical regions, which are identified 
by two natural conditions, to see clearly the temperature boundary. In the physical region 
the ratio appears a minimum value of 0.4 which is two times larger than that expected 
from the ideal behavior of the elliptic flow. The chemical potential effect, distinguished 
from the viscosity case, decreases the entropy density and thus contributes a final positive 
effect on the ratio. 

At last, we emphasize here that one should pay attention not only to the ratio of 
transport coefficients to the equilibrium entropy density but also to the ratio of transport 
coefficients to the non-equilibrium entropy density when to understand the transport 
or dissipative properties of QGP. Generally speaking, the entropy production dS comes 
either from the external sources d ex S and the internal dissipative and viscous processes, 
which is expressed as 

dS = d in S + d ex S. (54) 

As to the external sources, two basic aspects should be taken into account. One is the 
energy loss of parton when it is passing through the medium. The other might be the 
increase of the degree of freedom excited by the phase transition. Actually, these problems 
are really open and far from known, which need further study. 
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Appendix: Relaxation time 



In this appendix, we first define the relaxation time in Boltzmann equation, and then relate 
it with the shear viscosity which has been already obtained in the variational approach. 

The Boltzmann equation without time derivative and external force terms in the re- 
laxation time approximation is to substitute the collision term on the right hand side for 
the fluctuation of the distribution function scaled by the relaxation time r v 

t .BM*V>) = _SJ., (55) 
ax t v 

where /°(x;p) = [exp(— f3P v u u =p /i) ± l] 2 is the Jiittner distribution function [6]. With 
this definition and Eq. (jUJ), one could find 

S(T ZX ) = r v p j \^^\ {gfn f (p)[l - n f (p)} 

fdu z \ 

+9fnf(p) [1 - n/(p)] + Qb-ribip) [1 + n 6 (p)] } ( J (56) 

= (57) 

where the last line is the linear law. 

Noticing the velocity gradient is independent with momentum thus can be cancelled 
on both side, we can find out the relation between the relaxation time and the shear 
viscosity by trivially evaluating the momentum integral, 



3.246T 4 ( 1 + 0.246^ 



2 ' 

(5* 



21 



References 

[1] E.V. Shuryak and I. Zahed, Phys. Rev. D 70, 054507 (2004) 
[2] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 (2005) 



[3] B. Miiller, arXiv:nucl-th/0508062 and references therein. 



[4] P. Haensel, K.P Levenfish and D.G. Yakovlev, Astron. & Astophys., 381, 1080 (2002); 
ibid. 372, 130 (2001); ibid. 357, 1159 (2000) 

[5] Z. Xiaoping, K. Miao, L. Xuewen and Y. Shuhua, Phys. Rev. C 72, 025809 (2005) 

[6] S.R. Groot, W.A.van Leeuwen and Ch.G.van Weert, Relativistic kinetic theory 
(North-Holland Publishing Company, 1980) 

[7] A. Nakamura and S. Sakai, Phys. Rev. Lett. 94, 072305 (2005) 

G. Policastro, D.T. Son and A.O. Starinets, Phys. Rev. Lett. 87, 081601 (2001) 

[9] A. Buchel and J.T. Liu, Phys. Rev. Lett. 93, 090602 (2004) 

[10] PK. Kovtun, D.T. Son and A.O. Starinets, Phys. Rev. Lett. 94,111601 (2005) 



[11] K. Maeda, M. Nasuume and T. Okamura, ||hep-th /0602010 

[12] T. Hirano and M. Gyulassy, JHEP 0508, 086 (2005) 

[13] R. Hosoya and K. Kajantie, Nucl. Phy. B 250, 666 (1985) 

[14] S. Gavin, Nucl. Phys. A 435, 826 (1985) 

[15] P. Danielewicz and M. Gyulassy, Phys. Rev. D 31, 53 (1985) 

[16] D.W.von Oertzen, Phys. Lett. B 280, 103 (1992) 

[17] G. Baym, H. Monien, C.J. Pethick and D.G. Ravenhall, Phy. Rev. Lett. 64, 1867 
(1990) 

[18] E. H. Heiselberg, Phys. Rev. D 49, 4739 (1994). 

[19] P. Arnold, CD. Moore and L.G. Yaffe, JHEP 0011, 001 (2000) 

22 



[20] P. Arnold, G.D. Moore and L.G. Yaffe, JHEP 0305, 051 (2003) 



[21] Liu Hui, Hou Defu and Li Jiarong , arXiv: [hepjth/0508237, Eur. Phys. J. C 45, 459 
(2006) 

[22] S. Jeon, Phys. Rev. D 52, 3591 (1995). 

[23] E. Wang, U. Heinz and X. Zhang, Phys. Rev. D 53,5978 (1996) 

[24] E. Wang and U. Heinz, Phys. Lett. B 471, 208 (1999) 

[25] M.H. Thoma, Phy. Lett. B 269, 144 (1991) 

[26] H. Defu and L. Jiarong, Nucl. Phys. A 618, 371 (1997) 

[27] G. Aarts and J. M. Martinez Resco, JHEP 0211, 022(2002); ibid. 0402, 061 (2004); 
[28] G. Aarts and J. M. Martinez Resco, JHEP 0503, 074 (2005) 
[29] M. Basagoiti, Phys. Rev. D 66, 045005 (2002) 



[30] H. Defu, hep-ph/0501284 



[31] S. Jeon and L.G. Yaffe, Phys. Rev. D 53, 5799 (1996) 

[32] M.E. Carrington, H. Defu and R. Kobes, Phys. Rev. D 62, 025010 (2000); ibid. 64, 
025001 (2001) 



[33] Liu Hui, Hou Defu and Li Jiarong , arXiv: hep-ph/0602221 

[34] J.-P. Blaizot, E. Iancu and A. Rebhan, Phys. Rev. Lett. 85, 2906 (1999) 

[35] J.D. Brjoken, Phys. Rev. D 27, 140 (1983) 

[36] SART collaboration, Phys. Rev. Lett. 91,052303 (2003) 

[37] D. Teaney, Phys. Rev. C 68, 034913 (2003) 



23 



